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In today’s modern power devices and rapid growth power demands, the need 
for precise and accurate protection relays is a must for the power distribution 
system. That is, to segregate faulty sectors within fewer cycles, power relays 
should perform at the highest level of accuracy to detect abnormal conditions 
in power distribution. Therefore, this work will investigate the enhancement 
of the numerical relay testing in terms of harmonic distortions effect on the 
digitized output waveform as direct causes of relay failures. However, as it is 
an expensive process of testing the digitizing element of the numerical relay, 
this paper proposes a new algorithm of wavelet transforms in power quality 
signal processing testing using MATLAB simulation. As this newly proposed 
method of advanced waveform analysis algorithm will enhance the testing 
process of digitizing elements, and reduce data compiling complexity, a 
comparison between conventional fourier transforms testing and wavelet 
algorithm under abnormal conditions will be simulated based on inserting 
multi harmonics effect. As a result, based on the wavelet bank of filters, de- 


noising, and decomposition structure filters, wavelet has provided promising 
results in defining the effect of waveform distortion tripping time, fault 
location, total harmonic distortion, signal-to-noise ratio, and spurious-free 
dynamic range. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Emad Awada 

Department of Electrical Engineering, Al-Balqa Applied University Amman 
Amman, Jordan 

Email: emad.awada@bau.edu.jo 


1. INTRODUCTION 

In today's widespread population across the globe, and the fast revolution of technology from high 
voltage to low voltage power system networks, protection relay plays an important role in preserving the power 
system and eliminating abnormality within a few cycles [1], [2]. Yet, as power relay comes in different types and 
functions, such as current relay, and voltage relay, digital relays have developed more efficient functions in terms 
of current, voltage, and frequency measurement accuracy and power automation systems [3], [4]. Therefore, this 
work will investigate the numerical relay in terms of the internal logical microprocessor digitizing process as the 
bottleneck between the real and digitized spectrum to enhance the operation of power system protection. As 
known, the power system is subject to many interferences as an open surrounding power generation, transmission 
line, and distribution which result in a high rate of distortion fault occurrence [2], [5], [6]. In addition, with the 
significant growth of non-linear electronic types of equipment load, such as AC and DC variable speed drivers, 


Journal homepage: http://beei.org 


660 im) ISSN: 2302-9285 


switch mode power supplies, and power semiconductor controllers, harmonics are generated in voltage and 
current waveforms due to nonsinusoidal current draw and feed back into the power system. As a result, harmonic 
affects both electromechanical and digital relays by producing a pure signal beyond the fundamental frequency 
[7]-[11]. As a result, the North American electric reliability corporation (NERC) has mandated four reliability 
standards for the US government electronic return originator (REO) jurisdiction to insure power system protection 
and control. Standards such as; PRC-008 under frequency and Shedding, PRC-011 undervoltage load shedding, 
PRC-017 special protection, and testing program, and PRC-005 transmission and generation protection deal 
directly with relay operation and testing [12] as a major element in power system monitoring, measuring, and 
controlling protection process. Traditionally, static electro-mechanical relays were widely used to provide sole 
function characteristics [2], [3]. However, with the fast growth of power system generation and distributions, and 
dynamic changes in transmission and distribution grid, the technology of power relay technologies have evolved 
in the past few decades to provide multi-function characteristics. That is, by monitoring multiple input parameters, 
electronic relays based on digitized discrete electronic microprocessors were developed. Yet, as numerical relays 
algorithm structures are based on digitizing, and signal processing of analog current and voltage input waveform, 
relay digitizing operation constitutes a major data acquisition in terms of electrical quantities, accuracy, and speed 
[13], [14]. 

Many researchs have been conducted in the area of testing digital relay output waveforms to ensure the 
integrity of proper system performance. According to Wannous and Toman [3], a study was performed by 
investigating harmonics impact through the fast fourier transform (FFT) algorithm by injecting high values 
harmonics. Others had studied the practical effect of harmonic influences on electromechanical and 
microprocessor relays fundamental tone [7]—[ 10]. Some studies had explored the cost associated with quantifying 
the reliability and power quality as in [15]-[17] in terms of harmonics influences. However, the testing 
methodology of analog digitizers had shown a costly process in addition to some disadvantages as will in terms 
of a large number of samples, compiling duration, and calibration process [2], [3], [14], [18]. For instance, in 
testing digitizer characteristics, stimulus sine waveforms were primarily used as a simple mathematical model, 
and the easiest waveforms to generate frequencies with satisfactory reliability [19], [20]. Meanwhile, other types 
of waveform, square waves, for example, have been solely used for testing. As a result, in this work, the focus 
will be to investigate digital relay digitizer operation in terms of the real-time function waveform by implementing 
a new effective testing algorithm of wavelet transform into the microcontroller relay to ensure accurate 
measurement of power waveform dynamic operation and quantization process. 


2. THEORETICALCAL BACKGROUND 

In digitizing analog current or voltage waveform, several blocks of the filter process are used, to 
remove high-frequency components from the digitized waveform [3], [14], [19]. With a sampling frequency 
of 4,000 samples per second for 50 Hz frequency, IEC standard 61850, a large number of quantized data is 
inserted into the digital sinusoidal signal. Filters are adjusted to match amplitude angle using Fourier 
transformation as shown in Figure 1. Ideally, as the relay digitizer converts analog current and voltage 
waveform into the digital waveform, no additive or missing data due to quantization and sampling prosses are 
accounted for. 


AC Voltag e/ ial 

Current : Digital Numerical 
Analog —- Relay Digitizer Wave for Relay 
Input Microprocessor 


Figure 1. Digitizing process of the analog power waveform 


Yet, in a real-time application, generated digital signal would have an error embedded in the original 
waveform which affects system performance. As analog waveform X(t) converted to the digital domain, the 
output waveform X(n) can be presented with an error added as in (1). 


X (n) = X(t) +e (1) 
Where Xin(t) is the original signal and e is error value. 


Therefore, testing parameters of the relay digitizer should be a crucial factor to define the accuracy of 
the converted signal. Parameters such as multi-harmonics effect, signal to noise ratio (SNR), total harmonic 
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distortion (THD), and spurious-free dynamic range (SFDR) have a major effect on fundamental waveform [3], 
[18]—[20]. Noise and non-linear load distortion result in major harmonics alteration to fundamental waveform 
due to an increase in rms component value [3], [18], [21]. The power of these undesired digitized waveform 
components would affect the purity of the digitized waveform as in IEEE Std. 1241802.11g for SNR. That is, 
the ratio of signal power to the undesired signal components power, as noise and quantization error, as in (2). 


SNR(dB) = 20 log =2*“ms (2) 


Pnoiserms 


Meanwhile, as sensitive loads are affected by harmonic distortion, IEEE Std. 1241 defines THD as 
the power ratio of harmonics power summation to the fundamental waveform power as in (3). 


_ Pat P3t+P4t..+Poo 


THD (3) 


1 
For the current waveform, THD is the ratio of RMS summation to the fundamental waveform as in (4). 


: Yy=2 ems 
THD(i) = ~——™ « 100% (4) 


Irms 


Meanwhile, for voltage waveform components, THD can be shown as in (5). 
\2i=2 ve ms 
THD(v) = =e 8 100% (5) 


Where vj is RMS voltage. 

Such a phenomenon will degrade the original converted signal and lead to SFDR. That is, in IEEE 
Standard 1241, SFDR is the ratio of discrete fourier transform (DFT) regenerated digital signal amplitude to 
the largest harmonic (spurious) component observed over the full Nyquist band [19], [22], [23]. 


Amplitude of Fundamental (RMS) 
Amplitude of Largest Spur(RMS) 


SFDR = 20 log( ) (6) 


3. NUMERICAL RELAY OPERATION AND TESTING 

In the state-of-the-art power systems protection scheme, numerical relays are used with the advantage 
of programmable characteristics. Unlike conventional electromechanical and static relays, numerical relays do 
not require manual wiring and setting adjustments. Typically, numerical relay operation is based on measuring 
electric parameters in a network and converting measurements into digital mathematical and logical data to 
detect electrical faults and protect the network system. That is, a preset value of current or voltage needs to be 
defined as a threshold for numerical relay microprocessors to operate and disconnect from the electrical 
network [2]. Therefore, with quantizing binary inputs, configurable outputs, and programmable logic as shown 
in Figure 1, the digitizing process is an essential. Using a step-down voltage or current transformer, the 
secondary low side passes into a low-pass filter before digitizing which must have a high level of accuracy and 
sampling frequency at least twice the cycle of the highest frequency. Therefore, the proper digitization process 
of the numerical relay is very critical to ensur fundamental tone rather than distorting the waveform. 

As the conventional testing practice of analog digitizers based on waveform transformation, Fourier 
Transform is the core algorithm used in the frequency domain. Yet, as the fourier transform is based on the 
summation series of sine and cosine waveforms, characteristics of spectral impurity such as noise, harmonics, and 
frequency drifts may add and affect the digitizer testing output integrity [19]. As fourier transform articulates 
signals in the frequency domain, there will be no time representations in such an algorithm [19], [23]-[26]. 


4. METHOD 

In this work of testing the performance of numerical power relay digitizer, a model of relay harmonics 
effect was simulated based on distorted digital waveform. MATLAB simulation was used to generate a single 
phase power waveform, as shown in Figure 2, with the presence of noise distortion and hamanic effects. As 
power relay connected to power distribution system, four stages of convervation will be perform; i) 
voltage/current step-down transformation, ii) relay digitizer converter circuit (analog to digital), 111) digital 
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waveform data capture, and (iv) voltage/current data transformation for compiling process as shown in 
Figure 3. In this work, the simulated digitize output data was captured, analyzed, and validated through a model 
of conventional algorithm (fourier transform) and the new propesd algorithm (wavelet transform). Based on 
wavelet filter banks of denoising, decomposition, and localization in time and frequency domain, the behavior 
of power relay digitizer was monitored through extracting the detail components of output waveform and 
compiling SNR, THD, and SFDR as a major effect on numerical relay accuracy performance. 


Amplitude 


0 200 400 600 800 1000 1200 1400 1600 1800 2000 
Time 


Figure 2. Simulated relay digital output waveform 
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Figure 3. A newly proposed testing technique 


5. RESULT AND DISCUSSION 

As both algoritms of fourier and wavelet transformation were applied to analyze power relay output 
data, specific characteristics performance were tested, as follow, under abnormal and various fault types 
conditions. Several types of mother wavelet were applied to compare power relay digitizer fault mesearments 
with conventional algorithms in terms of SNR, THD, and SFDR. 


5.1. Signal to noise ratio 

In this simulation model, relay digital output waveforms were analyzed through FFT rather than DFT 
as a faster version. In addition, wavelet transform algorithm was used in computing SNR. By comparing the 
level of the signal components to the level of noise components in the signal background, the SNR of the 
digitized relay waveform was determined in terms of the ratio of signal power to noise power based on 
compiling all collected data samples in FFT as shown in Figure 4 and Table 1. 

Meanwhile, using the unique property of the wavelet decomposition process, providing significant 
information about frequency components, DWT was used to determine SNR based on much fewer samples. 
By decomposing the relay digitized waveform into multi-level and keeping the energy of the signal. In this 
work, 2" level decomposition was used to determine SNR based on detail coefficients cD2 as shown in 
Figures 5, 6, and Table 1. By comparing the computed SNR, it was clear that the Daubechie10 (Db10) had 
outperformed all other types of wavelet with SNR deviation of 1.35% from the FFT. 
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Figure 4. SNR of digitized relay waveform by FFT 


Table 1. SNR of digitized relay waveform by FFT and DWT 


Frequency FT Db12 Db10 Coefl Haar Bior3-1 
50 Hz 72.54 68.17 70.43 81.21 53.21 82.04 
60 Hz 71.48 67.89 69.75 94.31 51.82 86.95 
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Figure 5. SNR of digitized relay waveform by DWT 
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Figure 6. SNR reading illustration based on Table | 


5.2. Total harmonic distortion 


663 


As non-linear load feedback harmonics distortion into the power system’s, the perfect power waveform 
accumulat changes into the fundamental powert tone. Such accumulation leads to distorting power sinusoidal 
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waveform based on harmonic components magnitudes. In this work, both FFT as shown in Figure 7 and wavelet 
algorithms were applied to compute the effect of harmonics accumulation in terms of THD. 2™ level 
decomposition detail coefficients was used to determine the full scall range (FSR) of the relay output. By 
obtaining the summation square root of the differences between the detail coefficients and FSR, THD was 
measured based on the ratio of the square root to the maximum detail value (fundamental) as shown in 
Figures 8 and 9, Table 2. As a result, with 1.32% THD deviation from the FFT, Db10 had shown promissing 
results as comparing to other types of wavelet. 


4.5, THD = 82.35 


THD=85.58dB 
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Figure 7. THD of digitized relay waveform by FFT Figure 8. THD of digitized relay waveform by 
DWT 
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Figure 9. THD reading illustration based on Table 2 


Table 2. THD of digitized relay waveform by FFT and DWT 


Frequency FT Db20 Db10 Coef1 Haar Bior3-1 
50 Hz 85.58 82.35 83.06 88.30 87.65 88.92 
60 Hz 85.22 82.25 83.28 89.02 88.01 87.56 


5.3. Spurious-free dynamic range 

By defining the ratio of fundamental signal RMS to the worst spurious of interfered signals (noise and 
harmonics), SFDR distinguishes fundamental and largest component of unwanted interfering error. As a result, 
it is critical to determine SFDR in digital power relays to ensure proper operation response to grid frequency. 
In this simulation work, SFDR was measured based on 50 Hz and 60 Hz using both algorithms of FFT power 
spectrum and wavelet instantaneous dynamic range as shown in Figures 10 and 11 respectively. SFDR output 
results for FFT and wavelet analysis (2™ level detail coefficients), were illustrated in Table 3 and Figure 12. 
From Table 3 and Figure 12, once a gain the Db10 had outperformed all othe types of wavelet with SFDR 
deviation of 1.39% from the FFT. 
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Figure 10. SFDR of digitized relay waveform by Figure 11. SFDR of digitized relay waveform by 
FFT DWT 


Table 3. SFDR of digitized relay waveform by FFT and DWT 


Frequency FT Db12 Db10 Coefl Haar Bior3-1 
50Hz 12.27 11.37 11.96 15.44 14.46 15.98 
60Hz 13.23 11.92 12.85 16.23 16.08 15.96 
SFDR 
Dior) =r 
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Figure 12. SFDR reading illustration based on Table 3 


6. CONCLUSION 

As shown in this work, the digital power relay output waveform was tested to ensure proper 
performance of the digitizer and microprocessor for relays to protect transformers and power lines reliably 
when harmonics are involved. Both algorithms of FFT and discreet wavelet transform were used to compare 
testing results, performance time, data filtering, and compiling processes. Based on two fundamental 
frequencies, 50 Hz and 60, the sampling rate of 4,000 Hz, a total of 2,000 sample data points was collected and 
compiled based on FFT. Meanwhile, for the discreet wavelet decomposition, the 2" level of detail side had 
dropped the total number of data samples to 514 points. As a result, parameters of the relay output waveform 
such as SNR, THD, and SFDR were tested with fewer data samples to compile resulting in less complexity 
and shorter testing time, as will. In addition, based on these testing results, it was clear that daubechies wavelet 
(Db10), with average mesurment deviation of 1.35% from the FFT, was one of the most appropriate mother 
wavelet to conduct power waveform analysis based on the specific sinusoidal shape that match current and 
voltage sinusoidal waveform. 
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